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Abstract 

We consider a two-level system coupled to an environment that evolves non- 
adiabatically. We present a non-perturbative method for determining the per- 
sistence amplitude whose phase contains all the corrections to Berry's phase 
produced by the non-adiabatic motion of the environment. Specifically, it 
includes the effect of transitions between the two energy levels to all orders 
in the non-adiabatic coupling. The problem of determining all non-adiabatic 
corrections is reduced to solving an ordinary differential equation to which 
numerical methods should provide solutions in a variety of situations. We 
apply our method to a particular example that can be realized as a magnetic 
resonance experiment, thus raising the possibility of testing our results in the 
lab. 
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I. INTRODUCTION 



In the original Berry phase scenario fl[], the focus of attention is a quantum system with 
a discrete, non-degenerate energy spectrum. Its Hamiltonian H\R] is assumed to depend 
on a set of classical parameters R which represent an environmental degree of freedom to 
which the quantum system is coupled. The environment is assumed to evolve adiabatically. 
This produces an adiabatic time dependence in the quantum Hamiltonian, H = H\R(t)]. 
The time dependence of the quantum state \ip{t)) is determined by solving Schrodinger's 
equation using the quantum adiabatic theorem. Towards this end, one introduces the energy 
eigenstates of the instantaneous Hamiltonian H[R(t)}, 

H[R(t)]\E[R(t)}} = E[R(t)]\E[R(t)}} . 

It is further assumed that the environment is taken adiabatically around a loop in parameter 
space such that R(T) = R(0), and that the quantum system is initially prepared in an 
eigenstate |i?[R(0)]) of the initial Hamiltonian if[R(0)]. The quantum adiabatic theorem 
states that, at time t, the quantum system will be found in the state |.E[R(£)]) to within a 
phase factor, 



\ip(t)) = exp 



ij E (t) - \ f o dr£[R(r)]] \E[R{t))) . (1) 



The second term in the phase of the exponential is known as the dynamical phase and was 
already familiar from previous studies of the quantum adiabatic theorem. The first term 
represents Berry's discovery, and is referred to as Berry's phase, 

lEit) = i jf dr (E[R(t)]\^\E[R(t)}) . (2) 

In the cases where Berry's phase is physically relevant, 7^ is non-integrable: it cannot be 
written as a single- valued function of R over all of parameter space. Simon showed 
that the quantum adiabatic theorem has a line bundle structure inherent in it, and that 
Schrodinger's equation defines a parallel transport of the quantum state around the line 
bundle. Berry's phase is the signature that the associated connection has non-vanishing 
curvature. In this paper we will consider Berry's original scenario for a two-level system 
(2LS), though we will remove the adiabatic restriction on the environment. Our goal is to 
obtain the corrections to Berry's phase produced by non-adiabatic effects. 

The organization of this paper is as follows. In Section [IT], we introduce a non- 
perturbative method for determining all non-adiabatic corrections to Berry's phase. From 
the derivation, it will be clear that the effect of transitions between the two energy levels 
has been included to all orders in the non-adiabatic coupling. The problem of determining 
these corrections is reduced to solving an ordinary differential equation, to which numerical 



methods should provide solutions in a variety of situations. In Section [TTl| we work out a 
particular example in great detail. In Section |U1 A| we apply our method to this example and 
determine exactly the non-adiabatic corrections to Berry's phase. As a test of our method, 
in Section [111 B| we solve the Schrodinger equation exactly using a rotating frame transforma- 
tion, and use this solution to independently obtain the non-adiabatic corrections to Berry's 



phase. Comparison with the result obtained in Section HT A| shows that both methods yield 
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the same result. In Section [III C| we examine the corrections to Berry's phase obtained from 
our analysis in the limit of weak non-adiabaticity. We do this numerically and analytically, 
and compare our results with a previous result due to Berry. In Section |111D| , we discuss a 
magnetic resonance experiment that provides a realization of this particular example, and 
show how the non-adiabatic corrections to Berry's phase can be observed in measurements 
of the transverse magnetization. Finally, we make closing remarks in Section [TV. 



II. GENERAL ANALYSIS 

As mentioned in the Introduction, we consider a 2LS coupled to an environmental degree 
of freedom R(t) = R{t) (sin 9 cos (f>, sin 9 sin 0, cos^) with non-adiabatic time dependence. 
The coupling is described by the Hamiltonian, 

H(t) = R(t) a . (3) 

We denote the instantaneous eigenstates of H(t) by \E±(t)) with corresponding eigenvalues 
E±(t) = ±R(t). 

Because H(t) has non-adiabatic time dependence, transitions are possible between the 
two energy levels. Consequently, if we prepare the 2LS in the negative energy level \E-(0)) 
initially, there is a finite probability amplitude TL(t) to find the 2LS in the positive energy 
level \E + (t)) at time t. T_(i) is the transition amplitude, and the subscript indicates that the 
transition began in the negative energy level. Similarly, the probability amplitude that the 
system is found again in the negative energy level at time t defines the persistence amplitude 
P-(t). The subscript again indicates that the system was initially in the negative energy 
level. The amplitudes P+(t) and T + (t) have analogous definitions. 

The 2LS dynamics is determined by the propogator U(t,0) = exp[—(i/h) J * dr H(r)\: 

U(t,0)\E±(0)) = P±(t)\E ± (t)) + T ± (t)\E T (t)) . (4) 

We will determine P-(t) and T_(i) below, though our principle interest is in P-(t). The 
following derivation is easily adapted to determine P+(t) and T + (t), though we will not 
provide that derivation here. It proves convenient to write U(t, 0) as a 2 x 2 matrix: 

U(t,0)= £ U t ,\EM(EM = (T + + (t ) pj)) • 

To begin, we divide up the time interval (0, t) into n shorter time intervals of duration 
e = t/n by introducing intermediate times tk — he, (k — 0, • - • , n). Later we will let n — > oo. 
Clearly, U(k) = U(t}~, tk-i) propogates the state over the k-th sub-interval, and 

U(t,0) = U{n) ■ ■ ■ U{1) . (6) 

U(k) has the same structure as eqn. @: 

U(k)-( AP+{k) AT - {k) ) (7) 
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As the notation implies, AP±(k) and AT±(k) are the persistence and transition amplitudes 
corresponding to the k-th sub-interval. Noting that, 



le 



U(k)^l-"-H(k) + 0(e 2 ) 



(8) 



(E ± (k)\ ^ (E ± (k - 1)\ + e^(E ± (k - 1) 

and using eqns. (§), ©, and (§), one finds that, 

ie 



0{< 



AP ± (k) = l + ie^ ± (k)--E ± (k) 



(9) 



(10) 



; AT ± (k) = -eT ± (k) . 

Here a dot over a symbol indicates time differentiation, and 

ii±{k) = -(E ± {k)\E ± {k)) ; r ± (fc) = (£ T (fc)|£ ± (fc)) , (11) 

with 

r + (fc) = -r (fc) . 

7±(fc) are the Berry phases for the ± energy levels, and T±(k) are known as the non-adiabatic 
couplings for the ± energy levels. 

Inserting eqn. (^) repeatedly into eqn. @, and carrying out the necessary matrix multi- 
plications, one can show using induction that, 



P-(t) = l[AP4k) 



fe=l 



+ 



n Ap -( k ) 

k=n\+l 

+ ■ 



[AT + M] 



Til — 1 

n Ap +( k ) 

k=U2 + l 



[AT_(n 2 )] 



Vt 2 - 1 

n Ap -w 
. fc=i 



(12) 



We can make this equation more intelligible by introducing the amplitude V±(i,j — 1) to 
persist without any transitions in a given energy level over the interval {tj-i, t«): 



V±(iJ-l) 



UAP±(k) . 

k=j 



Eqn. (12) can then be re- written as, 

P_(t) = V-(n,0)+V-(n,n 1 )AT + (n 1 )V + (n 1 ,n 2 )AT_(n 2 )V-(n 2 ,0)+ ■■■ . (13) 

Each term in eqn. ( |T3"D gives the amplitude that the state of the 2LS follows a particular time 
sequence that begins and ends in the negative energy level. For any given time sequence, each 
sub-interval will have an amplitude associated with it which indicates whether a transition 
occurred during it (AT(k)), or not (AP(k)). Thus the first term gives the amplitude that the 
system undergoes zero transitions; the second term gives the amplitude that two transitions 
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occurred ( in sub-intervals n\ and 112). The remaining terms correspond to 4-transitions, 6- 
transitions, etc. . Only an even number of transitions are possible since the time development 
begins and ends in the negative energy level. Thus each transition out of this level must 
eventually be followed by a transition back into it. One can set-up a diagrammatic calculus to 
produce all the terms in eqn. (0), complete with rules for assigning a probability amplitude 
to each diagram, though we won't take the time to work that out here. 
Similarly, one can show that 

T_(t) =V+(n,m) AT_(ni)^-(7n,0) 

+V+(n,m) AT_(m) P_(m),7i2) AT + (n 2 ) V + (n 2 ,n 3 ) AT_(n 3 ) P_(n 3 ,0) 

+ ••• . (14) 

Here only an odd number of transitions can occur since the 2LS must finish in the positive 
energy level after having started in the negative energy level. Further simplification is 
possible if we introduce 

E(m,n 2 ) = AT + (m) V + (ni,n 2 )V-~ 1 (ni,n 2 ) AT_(n 2 ) , 
and insert 1 = V~(k)V-~ (k) appropriately into eqns. (fOp and (|14D. One finds that, 



and, 



P_(t) = V-(n,0) [1 + E(m,n 2 ) + E(n 1 ,n 2 )E(n 3 ,n 4 ) + 

2L(t) =7>_(n,0) {v+{n,n 1 )AT4n 1 )V-- 1 {n,n 1 )} 

x [1 + E(n 2 ,n 3 ) + E(n 2 ,n 3 )E(n 4 ,n 5 ) + 

Using eqn. (0), and recalling that n is large, one can show that, 
V±{n u n 2 ) = exp 



"1 



and, 



with, 



E(n 1 ,n 2 ) 



E €{iy ± {k)-{i/h)E ± {k)} 

k=n2 



[eF*{m)] [eF{n 2 ] 



F(m) = r_(m)exp 



T E c i(E + (k) - E.(k)) - h(j + (k) - 7-(*))} 



k=0 



So far, we have only considered one particular choice of intermediate times. We must now 
sum over all (maintaining the proper time orderings). This yields the following expression 
for F_(t): 



P_(t) = exp 



il-(t) 



h Jo 



(ItEJt) 



S(t) 



(15) 
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where, 



S(t) 



dyiF*( yi 



Mi 



+ 



rt ryi rx 1 ry 2 

I d Vl F*( Vl ) / dx 1 F(x 1 ) / dy 2 F*{y 2 ) / dx 2 F(x 2 ) 
o Jo Jo Jo 



and, 

F(t) = r_(t)exp 



dT5(r) 



E 4 



E-(r) 



h 



(7+(r) -7_(r)) 



(16) 



(17) 



We see that S(t) = A(t) exp[ip(t)] contains all the consequences of the non-adiabatic time 
dependence, and that it includes transitions between the levels to all orders in the non- 
adiabatic couplings r±(t). It is also clear that p(t) contains all the non-adiabatic corrections 
to Berry's phase J-(t). We close this Section by presenting a procedure for evaluating S(t) 
which promises to be useful in a variety of situations. 

It is a simple matter to write eqn. ([IB]) as an integral equation for S(t): 



S(t) = 1 



dyF*(y) / dx F(x)S(x) 



Introducing the auxiliary quantity I(t), 



I(t) 



dx F(x)S(x) 



S(t) 



1 dl 



(18) 



(19) 



F(t) dt ' 

and differentiating eqn. (|1^) with respect to t yields an ordinary differential equation for 
I(t): 



~d¥ 



di 

~dl 



+ \F\ 2 I = 



(20) 



the appropriate initial conditions are 1(0) = 0, and 1(0) = F(0) (note that 

either 



From eqn 

5(0) = 1 according to eqn. (|18 )). Determining I(t) reduces to solving eqn 
numerically or analytically. This is expected to be possible in a variety of situations. From 
I(t) we determine S(t), and hence, 

Im S(t) 



tan p(t) 



A(t) = ^(Re S(t)) 2 + (Im S(t) f 



(21) 



Re S(t) 

I(t) also allows us to express T_(i) more succinctly. Using the above results, one can show 



that, 



T-(t) 



exp 



i7+ (*) 



drEAr) 



I(t) . 



(22) 



Eqns. dl5|) — (p2|) constitute a general approach for determining completely the consequences 
of the non-adiabatic motion of the environment. Specifically, p(t) contains all non-adiabatic 
corrections to Berry's phase, while A(t) describes the reduced amplitude for the 2LS to be 
found in the negative energy level at time t due to transitions. In the following Section, 
we examine a particular example which is both experimentally realizable, and yet simple 
enough that our equations can be evaluated without approximation, and tested against the 
exact solution of the Schrodinger equation. 
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III. ANALYSIS OF A PARTICULAR EXAMPLE 



In this Section we will examine in great detail the interaction of a spin 1/2 with a time- 
varying magnetic field B(t). The magnetic field is assumed to precess about the z-axis 
at a fixed angle 9, at a constant precession rate <j)(t) = lu, and with constant magnitude 
|B(t)| = B. Such magnetic fields are encountered regularly in experiments involving nuclear 
magnetic resonance (NMR). As such, the results of this Section should be amenable to 
experimental test. The coupling of a spin to a magnetic field is described by the Zeeman 
Hamiltonian which, for a spin 1/2, has the same form as eqn. ([3D with the substitution 
R(t) = — 7^B(t)/2. Throughout this Section we will stick with the notation of eqn. ([$]), 
though it is a simple matter to substitute for R(t) when necessary. We will occasionally 
refer to R(t) as the magnetic field, though this is not literally true. 



A. Non-Adiabatic Effects: General Analysis 



In this subsection, we will determine all non-adiabatic corrections by evaluating S(t) 
using the general analysis of Section [Q|. To begin, we must determine the instantaneous 
eigenstates \E±(t)) of H(t). Straightforward analysis of eqn. (|3|) gives, 



\E + {t)) 



COS 2 



\E-tt)) 



sin 



(23) 



sin | exp [iut] I ' \ — cos | exp [iu>t] J ' 

where R(t) = i?(sin 9 cos ut, sm9smujt, cos#). Combining eqn. (|23|) with eqn. ( pTTj) gives 



F_(t) = -f sin0 



-iC 



F(t) = -iCexp[i5t] . 



\5(t) = f-ucos9 ^ 
F(t) is now inserted into eqn. ( p0|) to give: 

J - iSI + C 2 I = . 
This equation is easily solved since it has constant coefficients. One finds, 



(24) 



1(f) 



iuo sin 9 



■ exp 



i5t 



sin 



(25) 



with Vl = VS 2 + AC 2 . S(t) follows from eqn. (|T|): 



Re bit) = cos — cos h cos A9 sin — sin 

w 2 2 2 2 

lm bit) = — sm — cos h cos Ac^ cos — sin 

w 2 2 2 2 



(26) 
(27) 



Here cosA^ = 5/Qq: the physical significance of A9 will become clear in the following 
subsection. 
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B. Non-Adiabatic Effects: Rotating Frame Analysis 



In this subsection we obtain the exact solution to Schrodinger's equation for the particu- 
lar example considered in this Section, and derive from it S(t). This result will be compared 



with that obtained in Section |III A| , thus providing a test of our approach. 

The exact soultion can be found using a rotating coordinate frame analysis. In the lab 
frame, the Schrodinger equation is, 



(28) 



where H(t) is given by eqn. ([3D. We can transform to a frame that rotates with R(i) using 
the unitary operator, 



U(t) = exp 



iujt 

O'z 

2 



Writing \4>(t)) = U(t)\ip(t)), and substituting into eqn. 
in the rotating frame: 



gives the Schrodinger equation 



where, 



— , /{f2cos0-%4 

H = U^HU - mAu = 1 2 J 

Rsin9 



Rsin6 
{RcosO-^f} 



H is clearly time-independent, as expected, since the magnetic field is stationary in the ro- 
tating frame. The z-component of the magnetic field has been altered by the transformation. 
The magnetic field now makes an angle 9 with the z-axis given by, 



tan# 



iZsinfl 



RcosO - r -f 



(29) 



Note that = since H is real. 

Being time- independent, H has stationary states. The energies are 



E 



± 



± 



\ 



.Rcos 6 



^\ +^sin 2 fl 



hQ 



(30) 



and Q was defined in Sec. |III A|. The eigenstates are: 



\E+) 



cos • 



sin 



\E-) 



sm 



2. 

COS ' 



(31) 



The initial condition in the lab frame is | -0(0) ) = \E_(0)), and \E_(0)) is given in eqn. 
Since U(0) = 1, the initial condition in the rotating frame is ^(O)) = \E_(0)). Expanding 
|^(0)) in the basis \E±) gives, 



S 



1^(0)) = a + \E + ) +a_\E_) 
and application of the initial condition gives, 



A9 



A9 



sin 



cos 



(32) 



(33) 



Here A9 = 9 — 9, and is the same A9 as appeared in Sec. A| . One can see this by using 
eqn. (|^) and standard trigonometric identities to show that cos A9 = 5/Q , just as we found 
for A9 in Sec. [ill Aj . Physically, A9 is the change in the angle the magnetic field makes with 
the z-axis, as seen in the rotating and lab frames. From eqn. (BI3) we can immediately write, 



iQnt 



\E + ) + a_ exp 



iQnt 



\ip(t)) = a + exp 

Transforming back to the lab frame gives the exact solution \ip{t)) : 
\ip(t)} = exp — — 

r ;n .ii / cos f \ 

+ a_ exp 



E-) • 



x < a + exp 



2 



sin | exp [iut] 



2 



sm 



— cos I exp [iuit] 



(34) 



From eqn. (|3~4] ) we can obtain the persistence amplitude P-(t) = (E^(t)\ip(t)). Using 
eqn. (0), we find, 





iuot~ 




iQnt' 




'iQnt' 




exp 


. 2~_ 


|a^_ exp 


[ 2 


+ clL exp 


i 2 


} 



This expression can be straightforwardly re-written as, 



P_(i) = exp 



Z7_ 



iJl dTEJ - T) 



exp 



i5t 

T 



a, exp 



iQnt 



+ a_ exp 



iQnt 



Thus, the factor in curly brackets is S(t) (see eqn. (|i5|) ), as determined from the exact 
solution of Schrodinger's equation. The exact solution thus yields, 



He b (t) = cos — cos h cos A9 sm — sm 

w 2 2 2 2 



Im S(t) 



. St Q t , .„ 8t . Q t 

sm — cos h cos A9 cos — sm 

2 2 2 2 



(35) 
(36) 



Comparing eqns. (^) and ( |3"6D with eqns. ( PB] ) and P7|), we see that our approach gives 
precisely the same result for S(t) as the exact solution of the Schrodinger equation. 



C. Non-Adiabatic Effects: Numerical and Analytical Evaluation 

Here we explicitly evaluate the non-adiabatic corrections to Berry's phase. The exact 
result (including all non-adiabatic corrections) will be evaluated numerically. We also deter- 
mine an analytical approximation for these corrections in the limit of weak non-adiabaticity. 
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We shall see that the analytic approximation agrees quite well with the exact result in this 
limit. We also compare our results with an earlier result due to Berry. 

We begin by substituting eqns. fl2"6| ) and (|27| ) into eqn. (pT]). This yields, 



tan p 



#tan^ 



-tanf 



1 + 



#tan^ 



tan 



st 



(37) 



where g = cos = 6/Qq. Introducing e through the relation, 



et fl t 

tan — = q tan 

2 y 2 



(38) 



eqn. (^) becomes, 



tanp = tan 



(e - 5)t 



so that 



P 



(e-8) 



t + nil 



(39) 



We focus on the n = branch. Clearly, we must determine e. To do this, we differentiate 
eqn. (^) with respect to to at fixed t. This gives: 



de 
duo 



2 cos 2 f 



dg Q t gt dQ 
— tan 

duo 2 2cos 2 ^ duj 



(40) 



The appropriate initial condition for eqn. ([40]) is, 



2R 



s(uj = 0) = n (oo = 0) = — , 

h 

which follows from eqn. (|38|), since g(u = 0) = 1. It proves useful to write eqn. 
dimensionless form. For this purpose, we introduce the following definitions: 



in 



x 



2R 



T 



2R\ 



2R 



It also proves useful to define dimensionless versions of 5 and fl (see Sec. pi A|) : 

d = 5 = 1 — x cos 9 
e = (^) Q = Vl -2xcos^ + x 2 . 



(41) 



(42) 
(43) 

Thus, g = S/Q = d/e. We will integrate eqn. fl40| ) over the range < uj < ujf, and will use 
Tf = 2ir/u)f to define the time scale (t = (2n/ujf)s). Then, 

2ns 

T = . 

Xf 
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With all these definitions in place, eqn. ( fl0|) takes the following dimensionless form: 



de 
dx 



*f cos2 g 

SIT 



dg Tise 
tan h 



nsg 



de 



dx 



x f 



X f cos 



2 nse 



dx 



(44) 



Numerical integration of this equation yields e(x, t), which in turn gives p(t) via eqns. ( f41~|) 
and fl3~9"|). We remind the reader that this numerical result contains all non-adiabatic cor- 
rections to Berry's phase. The numerical result for p is given in Figure 1 (curve A). The 
integration was done with s — 1 (one precession cycle), 9 = 60°, and Xf — 0.3. 

It is possible to find a simple analytical expression for p(t) when x < 1. In this case, 
eqn. (f44|) can be treated iteratively. We only carry out the first step of this iteration pro- 
cedure since a numerical analysis is better suited to handle the case when x is not small. 
Noticing that e(0) = e(0), and since we are assuming that a; <C 1, the first iteration step 
substitutes e (x) = e(x) into the right-hand side (RHS) of eqn. (fill) . Solving the resulting 
differential equation yields ex, the improved e, which serves as the input for the second iter- 



ation step. Thus, e% is plugged into the RHS of eqn. ( |44l) which is solved to yield e 2 , etc.. 
Carrying out the first iteration yields the following differential equation for ex'. 

27TS 

dex _ de | x f sm ~ dg 
dx ^ dx 2tts dx 

Since S/ < 1, the sine function on the RHS oscillates rapidly, and so the second term on 
the RHS is not expected to contribute significantly to t\. We will see below that this is, in 
fact, the case. Thus we drop the second term and integrate the resulting equation to get, 

. . 1 f x f , de 
ei (x) = l + j / o dxg- . 



Since x <C 1, we can evaluate the integrand to second order in x using eqns. fl42l) and (^ 
and then carry out the integration. This gives, 

ei (x) = 1 — x cos 9 H sin 2 9 H sin 2 9 cos 9 + ■ • • 

u ) > 2 3 



From this we find, 



px (t) = cot cxx sin 2 9 + c 2 x 2 sin 2 9 cos 9 + 



(45) 



where Ci = 1/4, and c 2 = 1/3. Eqn. (|45|) is plotted in Figure 1 (curve B). We see that 
it agrees quite well with the numerical evaluation of the exact result (to within 1%) for 
the range of x-values considered. As anticipated, we see from Figure 1 that the discarded 
oscillatory term does not contribute significantly to p when a; C 1. 

Berry has worked out an adiabatic iteration procedure which generates an asymptotic 
expansion for corrections to Berry's phase . The procedure is adiabatic in that it ignores 
non-adiabatic transitions in all orders of iteration. This procedure yields a sequence of phase 
approximants, with each approximant containing powers of x to infinite order. Each itera- 
tion renormalizes the coefficients of the powers of x obtained in the previous iteration step. 
Berry has shown that the optimum number of terms to keep in the asymptotic expansion 
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is n ~ 1/x. The sequence of phase approximants initially improves with successive iter- 
ations, though ultimately, the sequence diverges because of transitions introduced by the 
non-adiabatic time dependence. 

Berry applied this procedure to the example we are considering in this Section. He 
carried out the first iteration step and worked out the resulting corrections to Berry's phase 
to order x 2 . His result has the same functional form as eqn. (^), though he finds c\ = 1/2, 
and C2 = 1. We believe that the discrepancy with our values for c\ and arises from 
truncation of the adiabatic iteration procedure at the first step. Such a truncation of the 
asymptotic expansion produces a non-optimum approximation for p(t) when i C 1 (see 
remark above), and consequently, non-optimum values for C\ and c<i- It is clear that Berry's 
intention was to illustrate his method; we believe that if a sufficient number of iteration 
steps were carried out, the two approaches would produce equivalent values for c\ and C2. 
The single-iteration result is plotted as curve C in Figure 1. 



D. Experimental Realization: Nuclear Magnetic Resonance 

One of the first observations of Berry's phase was by Suter et. al. using nuclear mag- 
netic resonance (NMR). In this experiment, the rotating magnetic field precessed about the 
z-axis in the manner assumed in this Section. Measurement of the transverse magnetization 
(M±(t)) allowed observation of Berry's phase. We now show that this same measurement 
(not so surprisingly) will also reveal the non-adiabatic corrections to Berry's phase deter- 
mined above. 

If initially the spin 1/2 has a component transverse to B(0), the spin will begin to precess 
about B(0). If B(t) does not evolve too rapidly, the spin precession simply follows B(i). To 
simplify the analysis, we assume, 

|V>(0)> = [|J? + (0)> + |£L(0)>] , 

corresponding to the spin being aligned initially along the x-axis in the lab frame. Using 
\ijj(t)) = U(t, 0)|^(0)), and eqn. (|), we have, 

\m) = j= [ {P+(t) + T_(f)} \E + (t)) + {P-(t) + T + (t)} \E4t)) ] . 

The transverse magnetization (M±(t)) = (M x (t) + iM y (t)) is given by 

(M ± (t)) = Tip d (t){ 1 M + } . (46) 

Here pd(t) = \ip(t))(ip(t)\ is the density matrix; I + = I x + il y is the raising operator for 
angular momentum; and 7 is the gyromagnetic ratio. We assume that t = lnjio so that 
\E±(2it/uj)) = \E±(0)). Also, in the basis \E±(0)), I + = |£+(0))(£_(0)|. Using these results 
in eqn. fl4"E| ) one finds, 

(M ± (2n/cu)) = y ( P - + T +) + T -*) 

= -L [P_P+* + T+P+* + P_T_* + T+T_*] . (47) 
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We have already evaluated P_. T_ can be evaluated using eqns. fl25|) and 
are determined by suitably adapting the analyses for P_ and TL. One finds: 



P_i_ and T_l 



P_ 
and 



A exp 



IT 



iuot 



1 + cos^) + ip 



Aexp 



iRt iuot 



T_ = —iC exp 



"~2~ 



ft 



1 — cos 9) — ip 



(48) 



(49) 



A is determined from eqn. and C from eqn. fl24j). 
Using these results in eqn. (PT7|) gives, 



(M ± (2ir/u)) 



I 11 ,2 

—A exp 



i8t + iu;t 



a; 2x 

— sin 2 9 H sin 2 9 cos 

2 3 



(50) 



We have assumed a; 1 so that we can use pi (i) from Sec [111 C to give an analytic approx- 
imation for p{t). In this limit, one can show that C 2 is of order x 2 , and so the second term 
on the RHS of eqn. (EDf) is negligible compared to the first. Thus, 



{M x (2im/uo)) 



-fhA 2 



cos 



2im 



1 



x 



cos 9 + x sin 9 



1 2 

— I — X cos # 

2 3 



(51) 



Here a; = u/jB when we substitute for R. The first two terms in the argument of the cosine 
function were observed by Suter et. al. [|]]. It would be interesting if this experiment could 
be repeated to look for the non-adiabatic corrections given in eqn. (|5~TD . 



IV. CLOSING REMARKS 

In this paper we have presented a non-perturbative method for determining all non- 
adiabatic corrections to Berry's phase. The problem of determining these corrections has 
been reduced to solving an ordinary differential equation (ODE) for which numerical meth- 
ods should provide solutions in a variety of situations. 

We applied our method to a particular example which can be realized as an NMR ex- 
periment, and whose Schrodinger equation can be solved exactly. For this example, our 
method could also be implemented exactly, and we saw that it yielded non-adiabatic cor- 
rections which were identical to those obtained from the exact solution. The non-adiabatic 
corrections to Berry's phase were evaluated numerically, and an analytical approximation 
was also obtained in the limit of weak non-adiabaticity. Both results were compared with 
a previous result by Berry. We also discussed how the non-adiabatic corrections to Berry's 
phase could be measured using NMR. 

We close with some final comments. (1) We stress that our method is non-perturbative. 
The object determined by the previously mentioned ODE contains non-adiabatic corrections 
to all orders in the non-adiabatic coupling. (2) The phase we determine is different from 
the Aharonov-Anandan phase B. In the scenario we consider, it is the system Hamiltonian 
which executes a cyclic evolution. Because the time dependence is non-adiabatic, the quan- 
tum system does not return to its initial state at the end of a cycle of the Hamiltonian, and 
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so its state will not, in general, execute a cyclic evolution. The phase we have evaluated 
is, in fact, the Pancharatnam phase 0-^], and we have explicitly seen that it reduces to 
Berry's phase in the limit where the non-adiabaticity goes to zero. We have also seen, for 
the example considered in Section |T|, that no geometric phase appears in the transition 
amplitude (see eqn. (pE5D), in agreement with Berry |J since <fi(t) is an odd function of t 
in this case. (3) It would be interesting to apply the method presented here to the case of 
an environment undergoing non-adiabatic stochastic motion. As discussed in Ref. [pj, the 
results of such an analysis should impact an ongoing controversy connected with the motion 
of vortices in superconductors. The controversy centers around whether certain Berry phase 
effects will be masked by a secondary process (connected with quasiparticle states bound 
to the vortex core) whose activation requires sufficiently large temperature and/or impurity 
concentration (see Ref. for further discussion and references). We hope to report on this 
application in a future paper. 
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FIGURE CAPTION 



Figure 1: A plot of the non-adiabatic corrections to Berry's phase p versus x = Tiuj/2R. 
Curve A is the numerical integration of eqn. (|44j ) which includes the effects of transitions 
between the two energy levels to all orders in the non-adiabatic coupling. In this calculation, 
s = 1, 8 = 60°, and Xf = 0.3. Curve B is our analytical approximation for p. Curve C is 
the first-iteration result of Ref. yj. 
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